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Physical properties of partially Ca substituted hole-doped Dy2Ti207 have been investigated by 
ac magnetic susceptibility Xac(T), dc magnetic susceptibility x{T), isothermal magnetization M{H) 
and heat capacity C'p(T) measurements on Dyi.8Cao.2Ti207. The spin-ice system Dy2Ti207 ex¬ 
hibits a spin-glass type freezing behavior near 16 K. Our frequency dependent Xac(r) data of 
Dyi,8Cao,2Ti207 show that the spin-freezing behavior is significantly influenced by Ca substitution. 

The effect of partial nonmagnetic Ca^"*" substitution for magnetic Dy®"*" is similar to the previous 
study on nonmagnetic isovalent Y®'*' substituted Dy2_a:Yj;Ti2 07 (for low levels of dilution), however 
the suppression of spin-freezing behavior is substantially stronger for Ca than Y. The Cole-Cole 
plot analysis reveals semicircular character and a single relaxation mode in Dyi.8Cao.2Ti207 as for 
Dy2Ti207. No noticeable change in the insulating behavior of Dy2Ti207 results from the holes 
produced by 10% Ca^"*" substitution for Dy^"*" ions. 

PACS numbers: 75.50.Lk, 75.40.Gb, 75.40.Cx 
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I. INTRODUCTION 

Rare earth pyrochlores R 2 S 2 O 7 {R = rare earth and 
B = transition metal) are known to present a wide range 
of interesting and unusual magnetic properties [ 1 - 11 ]. 
The three-dimensional network of corner-sharing tetra- 
hedra of R ions in these cubic symmetry (space group 
FdZm) pyrochlores cause highly frustrated magnetism 
making them a model material among the geometrically 
frustrated spin systems. An interesting situation is en¬ 
countered when the crystal electric field (CEF) forces the 
magnetic spin to point in a particular direction. The ob¬ 
servation of spin-ice behavior in Dy 2 Ti 207 and Ho 2 Ti 207 
is an exciting manifestation of such a CEF forced geo¬ 
metrically frustrated magnetism [3, 4, 7] . In these spin- 
ice systems CEF anisotropy forces the magnetic spins at 
the corners of tetrahedra to align parallel to the tetra¬ 
hedron axes (either towards the center or away from the 
center of tetrahedra). This causes an exotic magnetic 
state where even the ferromagnetic interaction between 
the Dy^+ (Ho^’*') magnetic moments becomes frustrated 
[3]. The energy minimum state of magnetic exchange 
and dipolar interactions is achieved by two spins point¬ 
ing inward and two spins outward on each tetrahedron, 
referred as the “2-in & 2-out” rule or “ice rule” [3]. The 
constraint condition on magnetic spins is analogous to 
the disordered protons in hexagonal water ice, accord¬ 
ingly the magnetic state of Dy 2 Ti 207 and Ho 2 Ti 207 is 
named as “spin-ice”. One of the major consequences of 
the ice rule is that a residual entropy remains at the low¬ 
est temperatures which is almost the same as the Pauling 
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value for water ice [4]. Another striking consequence is 
the realization of magnetic monopoles in spin-ice mate¬ 
rials which interact through the magnetic Coulomb in¬ 
teraction [ 8 ]. The diffuse neutron scattering study on 
Dy 2 Ti 2 07 shows evidence of Dirac strings and magnetic 
monopoles [9]. Polarized neutron scattering show evi¬ 
dence for a magnetic Coulomb phase in Ho 2 Ti 207 [11]. 
All these indicate a new state of magnetism in these ma¬ 
terials which led to intense research activities on R 2 B 2 O 7 
in recent years. 

Apart from the spin-ice behavior below T ^ 2 K, the 
ac susceptibility measurements on Dy 2 Ti 207 also show a 
frequency-dependent anomaly near Ty r; 16 K reflecting 
cooperative spin-freezing [12, 13], however the dc suscep¬ 
tibility or heat capacity shows no corresponding anomaly. 
This spin-freezing feature near 16 K is not seen in spin- 
ice Ho 2 Ti 207 which distinguishes the spin-dynamics in 
Dy 2 Ti 207 and Ho 2 Ti 207 . The spin-freezing behavior 
near Tf in Dy 2 Ti 207 is distinct from the spin-freezing 
in spin-glasses. In a spin-glass system the application of 
a magnetic held is found to quench the spin-freezing. In 
contrast, in Dy 2 Ti 207 the application of magnetic held 
tends to enhance the spin-freezing [ 12 ], thus indicating 
unusual dynamics of spin-freezing in Dy 2 Ti 207 . We fo¬ 
cus on this unusual glassiness in Dy 2 Ti 207 and study 
the ehect on spin-freezing behavior near Tf upon diluting 
the magnetic site and doping with holes, i.e. substitut¬ 
ing nonmagnetic Ca^’*' for magnetic Dy^+. The ehect of 
dilution on spin-freezing behavior near 16 K has been pre¬ 
viously investigated for nonmagnetic Y^"*" and Lu^“'" sub¬ 
stitutions [14, 15]. The investigations on Dy 2 _a;Ya;Ti 207 
by Snyder et al. [14, 15] revealed that the 16 K freez¬ 
ing is initially suppressed for low dilution concentrations 
{x < 0.4), however the freezing behavior was found to 
re-emerge for x > 0.4, clearly rehecting the inhuence of 
local environment on the freezing process. The 16 K 
freezing in Dy 2 Ti 207 is considered to be an Orbach-like 



2 


relaxation process which involves an excited CEF level 
for spin-lattice relaxations [16, 17]. The spin dynamics 
near Tf is controlled by transitions to/from the first ex¬ 
cited doublet. Therefore the effect of dilution at Dy site 
seems to result from a change in CEF. 

Here we report the physical properties of 
Dyi. 8 Cao, 2 Ti 207 based on ac magnetic susceptibil¬ 
ity Xac(T), dc magnetic susceptibility x(T), isothermal 
magnetization M{H) and heat capacity C'p(T) mea¬ 
surements. The results are compared with those of 
Dy 2 Ti 2 07 and the effect of partial nonmagnetic Ca^+ 
substitution for magnetic Dy^’*' on spin-freezing behavior 
near 16 K is addressed. The Xac(T) data show the effect 
to be similar to that of nonmagnetic substitution, 
however the suppression of spin-freezing behavior is 
much stronger in the case of Ca substitution. 


II. EXPERIMENTAL DETAILS 


Polycrystalline samples of Dy 2 _a;Caa;Ti 2 07 {x = 0, 
0.2) were synthesized by solid state reaction method. 
Dy 2 Ti 207 was synthesized using Dy 203 (99.9%, Sigma 
Aldrich) and Ti02 (99.995%, Alfa Aesar) in 1:2 molar ra¬ 
tio. The stoichiometric amounts were mixed and ground 
for half an hour using an agate mortar and pestle. The 
mixture was placed in a zirconium oxide crucible and 
fired thrice at 1350 °C in air with an intermediate grind¬ 
ing every 50 h, after which the finely ground mixture 
was pressed into a pellet and fired at 1350 °C for 50 h. 
For Dyi. 8 Cao. 2 Ti 207 synthesis Dy 203 , Ti02 and CaCOs 
(99.95%, Alfa Aesar) were taken in desired stoichiome¬ 
try, mixed and ground thoroughly, and fired at 950 °C 
for 50 h in air. The mixture was reground and fired 
at 1250 °C for 50 h with two more subsequent interme¬ 
diate grindings and firings at 1350 °C for 50 h each and 
then the sample was pelletized and fired again at 1350 °C 
for 50 h. A single phase sample of 10% Ca substituted 
Dyi. 8 Cao. 2 Ti 207 could be obtained this way. However, 
similar synthesis did not yield single phase sample for 
20% or 30% Ca substituted Dy 2 _xCaa;Ti 207 or 10% Na 
substituted Dyi. 8 Nao. 2 Ti 207 . A zirconium oxide crucible 
was used in each case. 

Crystal structure and quality of sample were deter¬ 
mined by powder x-ray diffraction (XRD). Bulk proper¬ 
ties were measured using facilities at Mag Lab, Helmholtz 
Zentrum Berlin. The dc magnetic measurements were 
carried out using a vibrating sample magnetometer 
superconducting quantum interference device (VSM- 
SQUID, Quantum Design) and the ac susceptibility was 
measured using the ACMS option of a physical prop¬ 
erties measurement system (PPMS, Quantum Design). 
Heat capacity was measured by the adiabatic relaxation 
method using the heat capacity option of PPMS. 
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FIG. 1. (Color online) Powder x-ray diffraction pattern of 
Dyi, 8 Cao. 2 Ti 207 recorded at room temperature. The solid 
line through the experimental points is the Rietveld refine¬ 
ment profile calculated for the Ca 2 Nb 207 -type face centered 
cubic (space group Fd3m) pyrochlore structure. The short 
vertical bars mark the Bragg peak positions. The lowermost 
curve represents the difference between the experimental and 
calculated intensities. 


TABLE I. Crystallographic and refinement parameters ob¬ 
tained from the structural Rietveld refinement of powder 
XRD data of Dy 2 Ti 207 and Dyi, 8 Cao, 2 Ti 207 . The atomic 
coordinates of Dy/Ca, Ti, Of and 02 atoms in space group 
FdZm are 16c (0,0,0), 16d (l/2,l/2,l/2), 48f (a:oi,l/8,l/8) 
and 8a (1/8,1/8,1/8), respectively. 



Dy2Ti207 

Dyi.8Cao.2Ti2 07 

Lattice parameters 

a (A) 

10.1231(1) 

10.1311(1) 

Hell (A^) 

1037.396(4) 

1039.843(6) 

Atomic coordinate 

2;oi 

0.4243(5) 

0.4254(5) 

Refinement quality 


1.86 

1.24 

Rp (%) 

2.11 

2.32 

Rwp (%) 

2.90 

2.96 


III. RESULTS AND DISCUSSION 

The powder XRD data were analysed by structural Ri¬ 
etveld refinement using the software FullProf [18]. The 
room temperature XRD pattern and Rietveld refinement 
profile for Dyi. 8 Cao. 2 Ti 207 are shown in figure 1. The 
almost single phase nature of the sample is evident from 
the refinement of XRD data. The structural refinement 
shows that no stuctural instability is caused by 10% Ca 
substitution for Dy, and like Dy 2 Ti 207 , Dyi 8 Cao, 2 Ti 207 
also has Ca 2 Nb 207 -type face centered cubic (space group 
FdZm) pyrochlore structure. The crystallographic pa- 
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FIG. 2. (Color online) The temperature T dependence of 
(a) real 'x and (b) imaginary x!' parts of the ac magnetic 
susceptibility Xac of Dyi,gCao. 2 Ti 207 measured at different 
frequencies from 100 Hz to 10 kHz in an applied ac magnetic 
field of 1.5 mT. 


rameters of Dyi. 8 Cao, 2 Ti 207 are listed in Table I and 
compared with those of Dy 2 Ti 207 . The crystallographic 
parameters are in good agreement with the values re¬ 
ported for Dy 2 Ti 207 [19]. 

The temperature T dependence of ac magnetic suscep¬ 
tibility Xac of Dyi, 8 Cao. 2 Ti 207 measured in an ac field of 
1.5 mT at different frequencies between 100 Hz to 10 kHz 
is shown in figure 2. It is seen from figure 2 that both 
the real part x' imaginary part x" of Xac show fre¬ 
quency dependent behavior. While at low frequencies 
there is no anomaly in x'{T) or x"{T), clear anomalies 
related to spin freezing are seen at higher frequencies. At 
10 kHz an anomaly in x'{T) is seen at around 19 K [fig¬ 
ure 2(a)] and that in x"{T) near 17 K [figure 2(b)]. This 
kind of frequency dependent behavior of Xac, i-O. ap¬ 
pearance of anomalies in x'i'^) x" at higher fre¬ 
quencies have also been observed in Dy 2 Ti 207 [12, 13]. 



FIG. 3. (Color online) The temperature T dependence of real 
x' and imaginary x" parts of the ac magnetic susceptibility 
Xac of Dyi,8Cao,2Ti207 and Dy2Ti207 measured at 10 kHz 
in an applied ac magnetic field of 1.5 mT. 


However, the frequency dependence is much weaker in 
Dyi, 8 Cao, 2 Ti 207 than that in Dy 2 Ti 207 . While at 1 kHz 
a well pronounced anomaly can be seen in the case of 
Dy 2 Ti 207 [12, 13], in the case of Dyi. 8 Cao. 2 Ti 207 no de¬ 
tectable anomaly is seen at 1 kHz in x'{T) and anomaly 
in x”{T) is very weak. The x'{T) and x''{T) data of 
Dyi, 8 Cao, 2 Ti 207 and Dy 2 Ti 207 at 10 kHz are compared 
in figure 3. Even though the anomaly temperatures of 
x'{T) and x"{T) are nearly same in both cases, the rel¬ 
ative changes in magnitudes at these temperatures are 
much weaker in Dyi, 8 Cao, 2 Ti 207 . The frequency de¬ 
pendence of Xa.c{T) for Dyi. 8 Cao. 2 Ti 207 is substantially 
weaker than the case of Y substituted Dy 2 _a;Ya;Ti 207 
[14]. Thus Ca^'*' substitution seems more effective in 
suppressing the spin-freezing behavior in Dy 2 Ti 207 than 
Y^+ substitution. 

The Cole-Cole plot of x" versus x' for Dyi. 8 Cao, 2 Ti 207 
is shown in figure 4 for T = 16 K. The Cole-Cole 
plot depicts the spin-dynamics at a given temperature 
and accesses the distribution of relaxation times of spin¬ 
freezing processes. While a semicircular character of plot 
indicates a single relaxation time or a narrow distribu¬ 
tion of relaxation times, the semicircle becomes flattened 
and/or distorted when there is a distribution of relax¬ 
ation times. From figure 4 we see that the Cole-Cole plot 
for Dyi. 8 Cao. 2 Ti 207 has a semicircular character similar 
to that of Dy 2 Ti 207 , though the center of the circle is not 
same for the two. Within the Cole-Cole formalism, the 
x"{x') data can be fitted to the equation of semicircular 
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FIG. 4. (Color online) Cole-Cole plot (imaginary part x” 
versus real part xO of fh® magnetic susceptibility for 
Dyi. 8 Cao. 2 Ti 207 and Dy 2 Ti 207 at different frequencies at 
16 K. 


arc given by [20] 

x"(x') 


Xo-Xs 


2 tan[(l —a)-|-] 

+ Xs)(xo - xO + 4tr n°;7i-5f] • 


( 1 ) 


Here xo is the isothermal susceptibility (in the limit of 
low frequency) and Xs the adiabatic susceptibility (in 
the limit of high frequency). The parameter a repre¬ 
sents the width of the distribution of relaxation times, 
a = 0 implies a single relaxation time and a = 1 an in¬ 
finitely wide distribution. The fits of x”(xO by ^0- (1) 
are shown by solid curves in figure 4. From the fits we 
obtained xo = 1.409(2) emu/mol, Xs = 0.49(3) emu/mol 
and (1 — 0 ) 77/2 = 89.46(2)° for Dyi, 8 Cao. 2 Ti 207 , and 
Xo = 1.552(4) emu/mol, Xs = 0.42(3) emu/mol and 
( 1 — 0 ) 77/2 = 89.46(3)° for Dy 2 Ti 207 . Accordingly we ob¬ 
tain a Ri 0.006 which is very close to 0 for both samples. 
The different values of xo and Xs reflect different center 
and radius of the semicircular arc for Dyi, 8 Cao. 2 Ti 207 
and Dy 2 Ti 207 . However, o is identical for both. The 
semicircular nature of Cole-Cole plot and the extremely 
small values of o suggest a single relaxation mode in 
Dyi. 8 Cao. 2 Ti 207 as well as in Dy 2 Ti 207 [12]. 

The zero-field-cooled dc magnetic susceptibility x)"^) 
of Dyi. 8 Cao. 2 Ti 207 measured \n H = 0.01 T and 0.05 T 
are shown in figure 5. The x)"^) data do not show any 
anomaly near Tf and paramagnetic behavior is evident 
from the xi'T) data at T > 2 K. The magnitude of x at 
low-T is quite large. A similar xi'^) behavior has been 
observed for Dy 2 Ti 207 [21]. From figure 5 it is evident 
that the x{T) measured at 0.05 T is higher in magni¬ 
tude than that at 0.01 T. This reflects the ferromagnetic 
nature of magnetic exchange as in the Dy 2 Ti 207 . At 
high temperature, the x(^) data of Dyi. 8 Cao. 2 Ti 207 fol¬ 
low Curie-Weiss behavior, xiT) = C/{T — Op). A linear 



FIG. 5. (Golor online) Zero-field-cooled magnetic suscepti¬ 
bility X of Dyi, 8 Cao. 2 Ti 207 as a function of temperature T 
for 2 < T < 25 K measured in magnetic fields H = 0.01 T 
and 0.05 T. Inset: Inverse magnetic susceptibility X~^(T) for 
2<T<50Kiniy = 1.0 T. The solid red line is the fit of 
the X~^{T) data by the Gurie-Weiss law in 20 K < T < 50 K 
and the dashed line is an extrapolation. 


fit of inverse magnetic susceptibility x~^{T) data (mea¬ 
sured m. H = 1.0 T) over 150 K < T < 300 K yields 
the Curie constant C = 28.22(5) emuK/mol and Weiss 
temperature Op = —14.9(6) K. The data when 

fitted over 20 K < T < 50 K yield Op = —0.55(1) K and 
C = 24.38(1) emuK/mol. The linear fit of the H = 1.0 T 
X~^(T) data in 20 K < T < 50 K is shown by solid red 
line in the inset of figure 5. Like Dyi. 8 Cao. 2 Ti 207 two 
different values of Op have also been found in the case 
of Dy 2 Ti 207 : dp = —15 K and —0.8 K from high-T and 
low-T fit ranges, respectively [22]. A positive dp has been 
reported for Dy 2 Ti 207 from the analysis of x{T) data be¬ 
low 40 K [4, 23]. Demagnetizing effect is argued to cause 
the dp to be negative and when corrected for this the 
dp would be positive such that effectively there is a weak 
ferromagnetic interaction in Dy 2 Ti 2 07 [21]. A similar sit¬ 
uation should be applicable to present data. Our sample 
for dc magnetic susceptibility measurements consisted of 
several small pieces of arbitrary shape and sizes therefore 
we do not have a good estimate of demagnetization fac¬ 
tor to correct the xi'T) data for the demagnetizing field. 
A crude approximation assuming spherical shape with 
demagnetization factor = 1/3 yields dp = 0.73(1) K 
on fitting the corrected x{T) data by Curie-Weiss law 
in 20 K < T < 50 K with C = 24.39(1) emuK/mol. 
Thus the demagnetization correction gives a positive dp 
and hence a ferromagnetic exchange. We estimate an 
effective moment fieB = 10.41 /ie/Dy from the value of 
C = 24.38 emuK/mol which is very close to the theoret¬ 
ically expected value of 10.61 fj,B for ®Gi 5/2 ground state 
of Dy^+ ions. 

The isothermal M{H) data of Dyi. 8 Cao. 2 Ti 207 mea- 
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FIG. 6. (Color online) Isothermal magnetization M of 
Dyi. 8 Cao. 2 Ti 207 as a function of applied magnetic field H 
for 0 < H < 7 T measured at the indicated temperatures. 


sured at selected temperatures of 2 K, 20 K, and 50 K 
are shown in figure 6. At 2 K the M{H) isotherm ex¬ 
hibits a rapid increase below 0.5 T and tends to satu¬ 
rate thereafter with a saturation value of magnetization 
Mg = 5.61 /re/Dy at 7 T. This value of Mg is close to 
the value expected for [100] direction Mg = 5.77/XB/Dy 
on account of FM interaction and Ising anisotropy [23]. 
However, the Mg value is much lower than the theoret¬ 
ical value of Msat = gjJ = 10 for free Dy^+ ions 
{gj = 4/3 and J = 15/2), in line with the presence of 
strong Ising anisotropy in this compound. The saturation 
behavior of M(H) isotherms depends on the tempera¬ 
ture: Mg decreases with increasing temperature. The 
overall M{H) behavior of Dyi gCao. 2 Ti 207 is similar to 
that of Dy 2 Ti 207 [21]. For a polycrystalline Dy 2 Ti 207 , 
Mg Ri 5^B/Dy at 1.8 K [15]. Thus we see an increase of 
about 12% in Mg value at 2 K with 10% Ca substitution 
for Dy. This effect of Ca substitution is quite distinct 
from Y substitution where dilution has negligible effect 
on saturation moment [15]. 

In Dy 2 Ti 207 there is a strong Ising anisotropy along 
the local (111) direction that determines the ice rule and 
hence the magnetic behavior of this compound. The 
M{H) data of Dy 2 Ti 207 obey ice rule (two-in, two- 
out conhguration) below 0.5 T and at fields higher than 
this Zeeman energy exceeds the effective nearest neigh¬ 
bor ferromagnetic interaction which causes a breaking of 
ice rule along the [111] direction [23]. The [111] direc¬ 
tion M{H) data for H > 2.5 T with Mg fv 5p,B/Dy 
at 1.8 K (much higher than the ice rule expected value 
of 3.33 ^u/Dy) reveal three-in, one-out spin configura¬ 
tion, clearly reflecting the breaking of ice rule at high 
fields [23]. For a spin ice system with two-in, two-out 
conhguration the predicted values of Mg for the [100], 
[110] and [111] directions are Msat/-\/3 (= 5.77^u/Dy), 


A7sat/V6 (= 4.08^B/Dy) and Msat/3 (= 3.33^B/Dy), 
respectively, where Mgat is the full saturation moment 
[23]. Therefore for two-in, two-out conhguration the 
weighted average Mg value for the polycrystalline sam¬ 
ple should be around 4.24/^B/Dy according to (Mg) = 
(6 M[iqo] -I-12 +8M[iii])/ 26, where the weight fac¬ 

tors 6, 12 and 8 are the number of equivalent directions 
for [100], [110] and [111], respectively. Our value of M 
for Dyi, 8 Cao, 2 Ti 207 is 4.04/rB/Dy at 0.5 T and 2 K 
which is compatible with the theoretically expected ice 
rule value [4.24 ^B/Dy] and is very close to the measured 
weighted average value (over [100], [110] and [111] di¬ 
rections) of M R! 4.0^B/Dy for single crystal Dy 2 Ti 207 
at 0.5 T and 1.8 K [23]. However, the observed M at 
0.5 T is about 14% higher than the M r 3.5p,B/Dy re¬ 
ported for polycrystalline Dy 2 Ti 207 at 0.5 T and 1.8 K 
in Ref. [15]. Thus, the saturation behavior of M{H) of 
Dyi. 8 Cao. 2 Ti 207 is substantially different from both un¬ 
doped as well as Y-doped Dy 2 Ti 207 [15]. While we do 
not have a clear explanation for the increased Mg for 
Dyi. 8 Cao, 2 Ti 207 we suspect that a change in crystal field 
effect or a change in exchange interaction by Ca substi¬ 
tution is the origin for the observed behavior. 

The charge/valence consideration would suggest three 
possible situations: 1) (Dy^+)i, 8 (Ca^+)o. 2 (Ti'*+) 2 ( 0^“)7 
i.e. with a deficiency of positive charge resulting in 
a hole doping. 2) (Dy3+)i,8(Ca^+)o.2(Ti^+)2(02“)7_5 
i.e. a dehciency of oxygen concentration to bal¬ 
ance the charge such that the composition is 
Dyi,8Cao,2Ti206.9 instead of Dyi. 8 Cao. 2 Ti 207 . 3) 
(Dy^+)i, 8 (Ca^''')o. 2 (Tii‘^+'^i''") 2 ( 0^“)7 i.e. a change in va¬ 
lence state of Ti, however this would require an increase 
in valence of Ti to compensate for loss in positive charge 
at Dy/Ca site if there is no change in oxygen concentra¬ 
tion. We would like to point out that Ti is not known to 
have a valence state higher than 4+, therefore this seems 
very unlikely. On the other hand from the view point of 
charge (valence) balancing an oxygen deficiency is pos¬ 
sible. However, an oxygen deficiency in Dy 2 Ti 207 has 
been found to reduce the saturation magnetization [24]. 
Contrary to such an expectation, an increase in satura¬ 
tion magnetization in Dyi, 8 Cao, 2 Ti 207 does not support 
the oxygen deficiency in this compound. An increase in 
saturation magnetization can be brought by oxygen defi¬ 
ciency if the deficiency concentration is more than what 
is required to balance the loss in positive charge upon Ca 
substitution such that a fraction of Ti becomes magnetic, 
i.e. Ti^"*" and the increased magnetization comes from Ti. 
Another possible origin of excess magnetization can be 
canting of Dy^+ moments from the easy (111) direction 
upon the partial removal of Dy^’*' spins from the vertices 
of a tetrahedron with two-in, two-out spin configuration 
as a result of Ca^+ substitution. Ca^+ is larger than 
Dy^+ in size and may distort the crystal field, changing 
slightly the Dy^"*" wavefunction and resulting anisotropy. 
If one of four Dy is replaced by Ca then such tetrahedra 
will have either two-in, one-out or one-in, two-out con¬ 
figuration. Thus assuming uniform distribution of Ca, 
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about 10% of the tetrahedra in Dyi. 8 Cao. 2 Ti 207 will no 
longer fulfil the two-in, two-out spin configuration and 
when these tetrahedra are subjected to high magnetic 
fields they may behave differently, canting of spins could 
be one such situation. 

The C'p(T) data of Dyi. 8 Cao. 2 Ti 207 are shown in fig¬ 
ure 7(a). It is seen from figure 7(a) that the C'p(r) data at 
T = 300 K is much lower than the classical Dulong-Petit 
value Cy = 3ni? = 33i? ~ 274.4 J/molK, where R is the 
molar gas constant and n = 11 is the number of atoms 
per formula unit [25, 26]. The Cp(T) data of Dy 2 Ti 207 
shown in figure 7(b) also fails to reach the high-T Dulong- 
Petit limit by 300 K. The lower value of Cp at T = 300 K 
than its high-T limit is consistent with the high Debye 
temperature (0d much greater than 300 K) of pyrochlore 
227 compounds. Like x(T), the C'p(T) data do not show 
any anomaly near Tf. At low-T an increase is observed in 
the Cp of Dyi. 8 Cao. 2 Ti 207 with decreasing T at T < 8 K 
[inset of figure 7(a)]. This behavior of C'p(T) is similar to 
that of Dy 2 Ti 207 [inset of figure 7(b)]. A comparison of 
Cp/T versus T plot for Dyi. 8 Cao. 2 Ti 207 and Dy 2 Ti 207 
in figure 7(c) shows that there is almost no change in 
upturn feature of C'p(T) by 10% Ca substitution. 

At low-T heat capacity is usually expressed by 
Cp{T) = yT -I- /3T^, where 7 and /3 are coefficients to 
the electronic and lattice contributions to heat capacity, 
respectively. The Sommerfeld coefficient 7 is nonzero 
for a metallic sample. However, for an insulating ma¬ 
terial 7 = 0, therefore C'p(T) can be described by j3T^ 
term only. We fitted the Cp(T) data of Dyi, 8 Cao. 2 Ti 207 
by C'p(T) = 7 T -I- /3T^ and found that 7 is indeed 
zero which is consistent with the insulating character¬ 
istic of Dyi. 8 Cao. 2 Ti 207 . The fit of Cp/T versus T^ 
by Cp/T = I3T^ in 12 K < T < 19 K is shown by 
solid red line in the inset of figure 7(c) which gives 
jS = 8.34(3) X 10“'^ J/molK"*. A similar analysis of 
C'p(T) of Dy 2 Ti 207 gives /3 = 9.78(5) x lO’^ J/molK^ 
in very good agreement with the reported value oi P = 
9.7 X 10“^ J/molK"* [27]. An estimate of Debye temper¬ 
ature 0D follows from /3, as 0d = (127r"‘ni?/5/3)^/^, ac¬ 
cordingly from the low-T C'p(T) data we get 0 d ~ 280 K 
for Dy 2 Ti 207 and 0 d ~ 295 K for Dyi. 8 Cao. 2 Ti 207 . 
The value of 0d for Dy 2 Ti 2 07 is in very good agreement 
with the reported value of 0 d = 283 K obtained from 
the C'p(T) data below 30 K [28]. 

The 0D values obtained above failed to fit the high-T 
C'p(T) data using the Debye model lattice heat capacity 
due to acoustic phonons therefore we analyzed the C'p(T) 
data using a combination of the Debye model and Ein¬ 
stein model of phonon heat capacity, thus accounting for 
both acoustic and optic modes of lattice vibration. The 
C'p(T) data were fitted by 

Cp{T) = mC'vDebye(T)-|-(l —m)C'vEinstein(T)-|-C'cEF(T) 

( 2 ) 

where Cy Debye (T) represents the Debye lattice heat ca¬ 
pacity, C'vEinstein(T) represents the Einstein lattice heat 
capacity and C'cef(T) represents the crystal electric field 
contribution to heat capacity. The phonon contributions 



FIG. 7. (Color online) Heat capacity Cp as a fnnction of 
temperatnre T for 1.8 K < T < 300 K measured in zero 
field for (a) Dyi, 8 Cao, 2 Ti 207 and (b) Dy 2 Ti 207 . The red 
solid curves are the sum of the contributions from the De- 
bye-|-Einstein models and the crystal electric field according 
to Eq. 2. Insets: Low-T Cp{T) for T < 20.0 K. (c) Cp/T 
versus T plot of low-T Cp{T) data of Dyi, 8 Cao. 2 Ti 207 and 
Dy 2 Ti 207 in 1.8 K < T < 20 K. Inset: Cp/T versus T^ plot 
for Dyi, 8 Cao. 2 Ti 207 for T < 20 K on an expanded y-scale. 
The solid red line is the fit to Cp/T = PT"^ in 12 K < T < 19 K 
and the dashed lines are extrapolations. 
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C'vDebye('r) and CyEinstein (T’) are given by [26] 

/ T \ ^ /•© d/t 

C'v Debye (r) = 9ni? — j (gX _ 

and 

/0p\^ g©E/T 

Cy Einstein(r) = 3 ni? 1 ^ —j (^^Ob/T _ X )2 > 

where 0 e is the Einstein temperature. For a two-level 
system C'cef( 7") is given by [26] 

Ccep(T) = r(^) (5) 


Since not much effect is seen on the Cp(T) data by 
10% Ca substitution for Dy, to analyse the Cp(T) data 
of Dyi. 8 Cao. 2 Ti 207 by Debye-|-Einstein model we used 
the same Ccef{T) that we estimated for Dy 2 Ti 207 . 
The sum of Ccef(T) and Debye-I-Einstein models fit for 
Dyi, 8 Cao. 2 Ti 207 is shown by the solid red curve in fig¬ 
ure 7(a). In this case we obtained m = 0.69(1), however 
the 0D = 766(9) K and 0 e = 172(4) K obtained from 
the fit are slightly different from those of Dy 2 Ti 207 . The 
increase in 0i3 and 0 e with Ca substitution is consistent 
with their formula mass dependence according to Eq. (6) 
as the formula mass is lowered by Ca substitution. 


IV. SUMMARY AND CONCLUSIONS 


where go is the degeneracy of the ground state, gi is the 
degeneracy of first excited state, and A is CEE splitting 
energy between the ground state and the first excited 
state. 

Under the action of CEE the (2 J -|- l)-fold degenerate 
ground state of Dy^+ {J = 15/2) splits into eight dou¬ 
blets, therefore go = gi = 2. There are few estimates of 
CEE splitting energy A for Dy 2 Ti 207 , however the CEE 
parameters have not been determined precisely by inelas¬ 
tic neutron scattering. The reported values of A are dif¬ 
ferent ranging from 100 cm“^ to 287 cm“^ [22, 29-31]. 
We use the value A = 100 cm“^ determined from the 
CEE analysis of magnetic susceptibility data [22]. The 
CEE contribution is only a few percent of the total Cp 
therefore uncertainty in estimate of Ccef( 7’) is not ex¬ 
pected to have significant effect on the analysis of the 
phonon contribution using Debye-|-Einstein models. 

The Ccef{T) was calculated using Eq. 5 for go = 
gi = 2 and A = 100 cm“^ and subtracted oft from 
the measured data and then the [C'p(r) — C'cef(U)] data 
of Dy 2 Ti 207 were fitted by Debye-|-Einstein models ac¬ 
cording to Eqs. (2)-(4). A reasonable fit is obtained for 
m = 0.71(1) i.e. with 71% weight to Debye term and 
29% to Einstein term. The sum of Ccef{T) and De¬ 
bye-I-Einstein models fit for Dy 2 Ti 207 is shown by solid 
red curve in figure 7(b). While fitting the data we used 
analytic Fade approximant for C'vDebye(U) as developed 
in Ref. [32]. The 0 d = 722(8) K and 0 e = 157(4) K 
were obtained from the fit. For Y 2 Ti 207 a value of 
0D = 967 K has been estimated [33] using the informa¬ 
tion from the first principles calculation of elastic prop¬ 
erties of Y 2 Ti 207 , which when scaled according to 


0D A 




1/3 

0DB 


( 6 ) 


gives 0D = 820 K for Dy 2 Ti 207 , where M and V are 
the mass and volume per formula unit, respectively and 
the subscripts A and B label the materials A and B, re¬ 
spectively. Thus the value of 0 d = 722(8) K determined 
from Debye-|-Einstein models fit for Dy 2 Ti 207 is consis¬ 
tent with the expected large value of scaled 0 d compared 
to the value of 0 d determined from coefficient /3 above. 


The effect of partial nonmagnetic Ca^"*" substitution 
for magnetic Dy^"*" in spin-ice system Dy 2 Ti 207 has 
been investigated by Xac{T), x{T), M{H) and C'p(T) 
measurements on Dyi, 8 Cao. 2 Ti 207 . The high-frequency 
Xac(T) data of Dyi.8Cao.2Ti207 show evidence for spin¬ 
freezing type anomaly near 16 K, however no correspond¬ 
ing anomaly is observed in x(T) and C'p(T) data. The 
x(T) and C'p(T) data of Dy 2 Ti 207 also do not show any 
anomaly near 16 K. The low-T Cp(T) data give Sommer- 
feld coefficient 7 = 0 which reflects an insulating ground 
state for Dyi. 8 Cao. 2 Ti 207 similar to Dy 2 Ti 207 . Thus 
the holes produced by 10% Ca^+ substitution for Dy^+ 
ions is not sufficient to change the electrical properties 
of insulating Dy 2 Ti 207 . The high-T C'p(r) data were 
analyzed by the Debye-|-Einstein models of lattice heat 
capacity along with a CEE contribution which yielded 
0D = 766(9) K and 0 e = 172(4) K for Dyi. 8 Cao, 2 Ti 207 
and 0D = 722(8) K and 0 e = 157(4) K for Dy 2 Ti 207 . 

The spin-glass type freezing behavior near 16 K in 
Dy 2 Ti 2 07 is strongly influenced by Ca substitution. The 
frequency dependent Xac(T) data of Dyi. 8 Cao, 2 Ti 207 
show that the suppression of spin-freezing behavior by 
Ca^+ is substantially increased compared to that by the 
nonmagnetic isovalent Y^+ in Dyi, 8 Yo. 2 Ti 207 [14]. An 
analysis of Xac{T) within the Cole-Cole formalism sug¬ 
gests that no change in semicircular nature of Cole-Cole 
plot of x"{x') results from Ca substitution and the spin 
dynamics of Dyi. 8 Cao. 2 Ti 207 is governed by a single re¬ 
laxation time like Dy 2 Ti 207 [12]. 

The Ca^"*" substitution for Dy from the viewpoint of 
charge (valence) balancing can result in a hole doping 
or an oxygen deficiency. Since an oxygen deficiency 
in Dy 2 Ti 207 has been found to reduce the saturation 
magnetization [24], the increased saturation magneti¬ 
zation in Dyi. 8 Cao. 2 Ti 207 is counter-intuitive for oxy¬ 
gen deficiency in this compound. As such the hole 
doping scenario seems more appropriate, however the 
holes produced by 10% Ca doping apparently fails to 
bring any noticeable change in band structure to reveal 
metallic nonzero 7 in heat capacity measurements. It 
should be noted though that we cannot completely rule 
out the possibility of O deficiency and the compound 



may have a composition Dyi, 8 Cao. 2 Ti 207_5 instead of 
Dyi.8Cao,2Ti207. 

In pure Dy 2 Ti 2 07 the spin relaxation has been sug¬ 
gested to exhibit an unusual double crossover, the spin 
relaxation is thermally activated at temperatures above 
a crossover temperature Tcross ~ 13 K below which the 
relaxation process is dominated by quantum tunneling 
down to Tice ~ 4 K and then again a thermally acti¬ 
vated behavior dominates below Tice as the spins become 
strongly correlated below this temperature [34]. Ehlers 
et al. [17] made an attempt to explain the difference 
between the 16 K freezing in Dy 2 Ti 207 (presence) and 
Ho 2 Ti 207 (absence) based on the competition between 
the thermally activated and quantum spin tunneling in 
these compounds. Considering that the spin dynamics 
near 16 K can be characterized by an Arrhenius type be¬ 
havior with an activation energy Ea ~ 220 K which is 
close to the CEF splitting energy (240-340 K), the spin 
freezing is considered to involve transitions to and from 
the first excited doublet [1, 17]. Since Dy^+ (J = 15/2) is 
a Kramers ion, a direct relaxation process is not allowed 
unless the degeneracy of ground state doublet is removed 
(e.g. by the application of magnetic field). Therefore an 
Orbach process is considered as an alternative mechanism 
of relaxation in Dy 2 Ti 207 which is a two stage process 
in which a spin first absorbs a phonon of the same mag¬ 
nitude as the energy of the excited crystal field state and 
then relaxes to a state lower in energy than its initial 
state by emitting a phonon of different energy than the 
absorbed one. Ho^+ (J = 8) on the other hand being 
non-Kramers ion, no Orbach-like relaxation is observed 
near 16 K in Ho 2 Ti 207 . Ehlers et al. suggested that 
Ho 2 Ti 207 can reveal 16 K feature under the application 
of magnetic field as confirmed by the neutron spin echo 
experiments [17, 35]. 

Considering the effect of Ca substitution on spin¬ 


freezing within the Orbach-like relaxation formalism that 
involves both the lattice degrees of freedom and the 
excited crystal field levels for spin-lattice relaxation, 
it seems that a significant change in the crystal field 
scheme/energy is brought by the Ca substitution. With 
a coordination number of 8, the (101.9 pm) has 
similar ionic radius as Dy^+ (102.7 pm), whereas Ca^’*' 
with an ionic radius of 112 pm is larger than Dy in size 
[36]. The larger size of Ca^+ might be at the origin for 
its effectiveness in changing the magnetic exchange in¬ 
teraction and/or crystal field, and could be the reason 
for stronger smearing of spin-freezing anomaly in ac sus¬ 
ceptibility than by Y^+. A comparative study of crys¬ 
tal field level schemes of pure Dy 2 Ti 207 as well as Ca- 
and Y-substituted ones would be informative in under¬ 
standing the magnetic behaviors of these compounds. In 
a recent study of effect of dilution on spin-ice behavior 
Scharffe et al. [37] found evidence for suppression of Paul¬ 
ings residual entropy in Dy 2 _a;Ya;Ti 207 for x > 0.4 with 
a non-degenerate ground state and accordingly a rapid 
crossover from spin-ice behavior in Dy 2 Ti 207 to a weakly 
interacting single-ion physics is inferred for intermediate 
dilution in Dy 2 _a;Ya;Ti 207 . It would be of interest to see 
how Ca substitution for Dy influences the spin-ice behav¬ 
ior and associated monopole dynamics. Further investi¬ 
gations in this direction would be enlightening. 
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